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Abstract 

In the cluster expansion framework of Bose liquids we calculate ana- 
lytical expressions of the two-body, three-body and four-body diagrams 
contributing to the g.s. energy of an infinite system of neutral alpha- 
particles at zero-temperature, interacting via the strong nuclear forces 
exclusively. This is analytically tractable by assuming a density depen- 
dent two-body correlation function of Gaussian type. For the a — a po- 
tential we adopt the phenomenological Ali-Bodmer interaction and semi- 
microscopic potentials obtained from the Gogny force parametrizations. 
We show that under such assumptions we achieve a rapid convergence in 
the cluster expansion, the four-body contributions to the energy being 
smaller than the two-body and three-body contributions by at least an 
order of magnitude. 

Key words: Equation of state, nuclear matter, supernova explosion, clus- 
ter expansion, quantum liquids. 

1 INTRODUCTION 

Understanding the properties of a matter has retained a lot of attention 
in recent years. This situation is mainly due to the believe that this type of 
hadronic matter occurs in astrophysical environment in deconfined form. In 
the debris of a supernova explosion, a substantial fraction of hot and dense 
matter resides in a particles and therefore the equation of state of matter 



1 



at subnuclear densities is essential in simulating the supernova collapse and 
explosions and is also important for the formation of the supernova neutrino 
signal pQ. 

The aim of the present work is to investigate the equation of state of a 
matter from the standpoint of the cluster expansion method of Bose liquids. 
We consider a cold (T=0) system of a-particles interacting only by means 
of the strong nuclear force. Similarly to the case of ordinary nuclear matter 
(composed of protons and neutrons) the Coulomb interaction is switched-off. 
The internal structure of the a clusters is accounted only in the determination 
of the a — a potential, the single particle structure being incorporated in the 
cluster densities that are folded with the effective nucleon-nucleon (NN) inter- 
action. In continuation to the previous assumption, no Pauli blocking effects 
are included. Naturally, since the constituents of the a particles are fermions, 
one should expect that the Pauli principle is manifest when two or more a 
clusters start to overlap. As revealed by the work of Ropke and collab. [2] one 
should expect from the action of this principle a dissolution of the a cluster 
in protons and neutrons above the so-called Mott density, which presumably 
lays between a fifth and a third of the nuclear matter saturation density. It 
is well known !3J. that two-body correlation functions (TBCFN) obtained by 
minimization of the energy functional truncated at the lowest order may lead 
to an unphysical deep minimum. We adopt the cluster expansion method of a 
Bose liquid [I] and use the simple Jastrow ansatz involving state-independent 
two-body correlation functions . These are taken in a Gaussian form, without 
overshooting near the healing distance, 



The parameter (3 is determined from the normalization condition for the cor- 
relation function [5j 



This condition ensures that the mean square deviation of the correlation func- 
tion from unity is a small quantity and has an exponential healing. Con- 
sequently the dependence j3 = f3(p) , where p = N a /V is the density of 
a-particles, reads, 



For the a — a interaction, we adopt two types of Gaussian-like potentials 
containing a short range soft repulsive part and a long range shallow attractive 
part. The first one is the S-state Ali-Bodmer (AB) potential [6]: 








I 1/3 



(3) 




(4) 
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In this expression, V A = 130 MeV, V R = 475 MeV, \i A =0.475 fin" 1 , p R =0.7 
fm _1 . This potential obtained by a fit of the low energy a — a phase shifts, 
can be considered as an approximation to the super symmetric partner of the 
deep potential of Buck et al. [7]. The second type of potential is a sum of 
three Gaussian and is derived from two recent parametrizations of the Gogny 
effective N — N force [8]. Two explicit forms, labeled (Dl) and (DIN), are 
given in [5]. We have checked that these potentials satisfying the integrability 
condition / °° r\v(r)\dr < oo are in agreement with the Levinson theorem in 
the sense given in [TU]. In the absence of the Coulomb interaction, the a — a 
interactions (AB) and (DIN) provide a weakly bound g.s. J w = + . In 
contrast to the aforementioned interactions, the a — a (Dl) interaction is 
characterized by the absence of bound states. 

In section 2 we apply the cluster expansion method to our problem and 
express all terms of the development in a compact form. Our results are 
discussed in section 3. 



2 CLUSTER EXPANSION METHOD 

Let the energy of a system of strongly interacting iV-bosons in the Jackson- 
Feenberg form |llj . 

E= -pN J drg(r)v*(r) , (5) 

where g(r) is the radial distribution function and v* is the effective Jackson- 
Feenberg potential, 

h 2 

v*(r)=v(r) V 2 ln/(r) . (6) 

2m a 

In the case of a-matter the energy Eq. ([5]) is measured relative to the rest 
energy of a free a-particle. Using Eqs. (Q]) and (UD we obtain 



2c/3 2 (3-e^ 2 (3-2r 2 /3 2 )) 



„2, ,2 „ _^2,,2 



v*(r) = K — —2 '- + V R e~ r »n - V A e~ r »a , (7) 

(1 - e r P ) 

where c = h 2 /2m a . The main trick allowing the cluster expansion of the above 
expression consists in expanding the radial distribution function in powers of 
the small parameter to = J drh{r). It uses the fact that the function h = f 2 — 1 
is of short-range. Accordingly, the cluster expansion of g(r) reads, 
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g(ri 2 ) = f 2 (r 12 ) |l + p J dr 3 h(r 13 )h(r 23 ) 



h(ri 3 )h(r M )h(r u ) + 2h(ri 3 )h(r 24 )h(r u )h(r u ) 



+ 0(p 3 )\ . (8) 



+ p 2 J dr 3 J dr 4 

+ ]jh{ri 3 )h{r 23 )h(ri 4 )h(r 24 ) 

+ ^K r n)h{ru)h{r 23 )h{r 2 A)h{r 3 A) 

A diagrammatic expansion of the radial distribution function is depicted in 
Fig. [TJ For a given n-body diagram there are field points (open circles) and 
dummy points (filled circles). For each dummy point there is an integration 
p J dfi. A bond between open points involves a factor f 2 in the integrand. Any 
other bond implies a factor h = f 2 — 1 in the integrand. Mutatis mutandis, the 
ground state energy per a-particle is then expanded in powers of the density : 

E = E 2 + E 3 + E 4 + ... (9) 

Above E n stands for the contribution to the energy per particle arising from 
n-body diagrams. In what follows we calculate these first three terms of the 
expansion. 

Figure 1: Diagrammatic expansion of the radial distribution function. There 
are 4 independent diagrams contributing to the 4th order term: ring, diagonal, 
opened and connected. 

Partial contributions to the total energy are displayed using the a — a 
Ali-Bodmer interaction (Eq.(U|)). 

2.1 Two-body diagram E 2 

Let us split E 2 only into kinetic and potential components : 
IP 



e 2K = 2nJ dr ^ dr 2f 2 ( r i2) 



2m n 



V^ln/(n 2 ) 



(10) 



E 2 y = dr i dr 2f(n2)v(r 12 ) , (11) 

where f2 is the integration volume. Consider the general 6-dimensional integral 

h = ~ J d ri dr 2 p(r 12 ) , (12) 
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where r\2 = r\ — r<i and p a generic function. With the unitary transformation 
(unit Jacobian), 

r 1= R+ls, (13) 



To = R S 

2 



the integral I2 (Eq. (fT2|) ) reads: 



1 p f r°° 
I2 = - — / dRds p(s) = 2np / ds s 2 p(s) 

2 u J Jo 



(14) 



(15) 



Note that integration over the cm. variable R gives the integration volume 
fi. Therefore we have 



E 2 K = -\pjds f 2 (s) 



2m, 



-V 2 ln/(s) 



2m r 



f(s)V 2 f(s) - (Vf(s)Y 



E-2V = \pj ds f(sMs) ■ 
Summing up these two contributions we obtain, 



E2 = ^P I drf z (r)v*(r) 



or in analytical form 



1^/2 



/ 3\/2c 



+ 2V R F R - 2V A F A 



1 



Defining the auxiliary function 

Xij(P,p) 
we then have 

Fi = -Xo,i{/3,Pi) + 2xi,i(P, Hi) - X2,i(P,p-i) 



(i/3 2 + jh 



2^§ ' 



(16) 
(17) 

(18) 

(19) 

(20) 
(21) 



where (i = A,R). The dependence of E2 on density is displayed in Fig. 
[2j We observe that already the component E2 has a shallow minimum at 
a density almost two times the saturation density of normal nuclear matter. 
This is in contrast to the result of ref. [9], where it has been shown that 
TBCFN's obtained from Pandharipande-Bethe equation lead to a collapse of 
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Figure 2: Density dependence of the E 2 component. 



E 2 component. This effect arises entirely from the density dependence of our 
particular functional form of TBCFN and its derivatives. 



2.2 Three-body diagram E 3 



The diagram corresponding to the three-body energy is given on the left 
panel of Fig. [3l By definition, 



ti 2 



2m a 



V?ln/(n 2 ) 



Ezv = \ (j- ) j dr 1 dr 2 drzf 2 {r ) > )h(r- 2 :\ )/'( <:!i )v{r [2 ) 



(22) 



(23) 



Like previously for the E2 component we consider the following generic integral 
h = ^ (jA J dr 1 dr 2 dr 3 p 1 (r 12 )p2{r23)P3(r3i) 



Figure 3: E 3 diagram and <5-diagram 



1 (P r 



n 



l[ dn J] dR iPl (R 1 )p2(R 2 )P3(R3) 
i i 

x5(Ri - n + r 2 )5(R 2 - r 2 + r 3 )5(R 3 -r 3 + n) . (24) 



Integration of variables is straightforward and we obtain, 

i*3 = \p 2 J dR 1 dR 2 dR 3Pl (R 1 )p 2 (R 2 )p 3 (R 3 )5(R 1 + R 2 + R 3 ) . (25) 

Eq. (|25|) is useful for both numerical and analytical integration since the angu- 
lar dependence is isolated in a 5 function as shown in the right panel of Fig. 
[3l Introducing the double-folding (5-integral, 

Vl q {R x ) = J dR 2 dR 3 p(R 2 )q(R 3 )5(Ri + R 2 + R 3 ) , (26) 

the integral I 3 is reduced to, 

I 3 = l -p 2 J dR lPl (R 1 )V^ p3 (R 1 ) . (27) 

Applying this technique to the E 3 term, we have, 

E 3 k = \ P 2 J dR 1 f(R 1 ) [f(Ri)V 2 f(Ri) - (Vf(Ri)) 2 } vUr,) 

(28) 

£sv = \p 2 J dRtfiRMRlWlhW ■ ( 29 ) 
Consequently the total three-body energy contribution reads, 

E 3 = E 3K + E 3Y = X -p 2 J dR 1 v*{R l )f 2 {R 1 )Vl h {R x ) , (30) 

where we have defined according to the prescription (|26p 



Vl h (Ri) = J dR 2 dR 3 h(R 2 )h{R 3 )5(R l + R 2 + R 3 ) . (31) 



For our particular selection of the TBCFN, 



TT 



3/2 / 1 



-H 2 



4^/3 



-e s 



+ V2, 



e 2 



and 



7i 3 p 2 r 



0» 



wi J + w * " VaGa . 



(32) 



(33) 



where for (i = A,R), 



Gi = —[xiAPiVi) ~ 2 X2,i(/3,Mi) + X3,xW,fM) + 32(xi,2(/3,^i) - 2x 3 , 2 (/3,Mi 



+ X5,2(P> Mi) - X2,3(/3, Mi) + 2X5,3(^, Mi) ~ X8,3(/?, Mi))]- 



(34) 




Figure 4: Density dependence of the £3 component. 



The density dependence of the three-body energy is given in Fig. HJ We 
notice that contrary to E2 which is attractive, £3 is only weakly attractive at 
low density and becomes strongly repulsive with increasing density. 
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Figure 5: (A) Ring, (B) diagonal, (C) opened and (D) connected diagrams 
contributing to the E 4 component. 

2.3 The four-body diagram E 4 

The diagrams contributing to the E4 component are depicted in Fig. [5j 
They are dubbed as ring, diagonal, opened and connected diagrams. Explicit 
expressions of these components are, 

1 o 3 f 4 

Ear =2^7 U d ^ v *(n2)f 2 (r 12 )h(r 23 )h(r 34 )h(r 41 ) , (35) 
1 o 3 f 4 

E m = 2-^ J ]Jdr i v*(r 12 )f 2 (r 12 )h(r 23 )h(r 34 )h(r 41 )h(r 24 ) , (36) 

E 4 o = \~ J Yldnv^r^f^ruMr^hirMr^hirrs) , (37) 

E 4C = \~ J II^^*^)/ 2 ^)^^)^^)/!^!)^^)^^) • (38) 



2.3.1 Calculation of E 4 r ring diagram 

It is easy to see that the ring diagram contribution to the energy is given 
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Figure 6: Density dependence of the -E4R component for the AB potential. 

by, 

1 r 4 

E 4R = -p 3 J YldR l f 2 (R 1 )v*(R 1 )h(R 1 )h(R 2 )h(R 3 )h(R 4 )5(R 1 +R 2 +R3+R4) 

(39) 

or equivalently, 

2 

Ear = \p 3 JfldR i f 2 (R 1 )v*(R 1 )h(R 2 )Vl h (R 1 + R 2 ) 

3 

= 2 p3 J Tl dR *f 2 ( R i) v *( R iMR2)Vl h (R 3 )8(R 3 -Rt- J?a),(40) 

in terms of the function V£ h defined in Eq. (|31|) . In the hierarchy of 5 kernel 
let us introduce, 

V p S *, r (Ri) = J dR 2 dR 3 p(R 2 )V q 5 r (R 3 )5(R 3 - R 1 - R 2 ) . (41) 
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The final integral is simply, 

E,r = \ P 3 J dR 1 v*{R 1 )f\R l )Vl% >h {R l 



(42) 



Using the expression for V£ h given in Eq. ([32|) we obtain the expression of 



V 5S 



7T 3 / 1 



/? 6 V 24^/3 



e 3 



12 -2bV 



R 2 p 2 
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5^/5 3^/3 
The corresponding expression for is: 



8^/2 

-R 2 I3 2 

e s 



(43) 



vr 9 /y 



C /3 3 
+ 



+ 
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In the equation 
1 



P \4 128^/2 4^/3 
240 \ 

the iJj's are defined by 



+ 



(44) 



Hi 



16 



(-6xi,2 (A Mi) + 12X3,2 03, w) - 6x5,2 (/?, Mi) - 64xi,3 (A Mi)+ 



X2, 3 (/3,Mi) + 128X4,3(0, Mi) " 2X5, 3 (/3,Mi) " 64 X 7, 3 (/3, Mi) 
+X8,3(/3, Mi) + 96x2,5(/3, Mi) - 192x7,5 (P, Mi) + 96xi2,s(/3, Mi)) , 

where (i = A,R). 

From the inspection of Figj6] we infer that the ring contribution to the 
four-body energy has the same behavior as E% save for a factor -1. 



2.3.2 Calculation of E^q diagonal diagram 

In order to isolate the angular variables, we proceed as follows, 

1 3 r 4 

E 4 D = ]Jdr i v*(r 12 )f 2 (r 12 )h(r 23 )h(r u )h(r 41 )h(r 24 ) 

E 4 d = 2~£ f\{dR i Jldr i v*{R l )f{R{)h{R 2 )h{R^h{R A )h{R h ) 
1 u J l l 
x5(Ri —r\ + r 2 )5(R 2 - r 2 + r 3 )5(R 3 - r 3 + r 4 ) 

x5(R4 -r A + r 1 )5(R 5 -r 2 + r 4 ) . (45) 
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Performing the same manipulations as above we obtain, 
1 r 5 

^4D = 2 • -p 3 y l[dR i f 2 {R 1 )v*{R 1 )h(R 2 )h{R 3 )h(R 4 )h(R 5 ) 

8(R 5 + R X + R 4 )5(R 5 -R 2 - R 3 ) ■ (46) 

We observe that variables Ri,i = 1,2,3,4 are decoupled. It can be seen 
easily that: 

£ 4D = P 3 J dR 1 dR 4 dR 5 f 2 (R 1 )v*(R 1 )h(R 4 )h(R 5 )Vl h (R 5 )S(R 5 + #i + R 4 ) . 

(47) 

Finally we have, 

7T 9 / V (c (27 26624^/3 - 7344 2432^6 - 10368 
4D " 288^ U V ¥ 49^7 + 25V5 

945 mi r 15552 580608 442368 \ T _ "\ 

I 74^ - — j= + — ._ - — ._ + Vkl>fl - Va Ai , (48) 



2^ * 121v/TT 169\/13 361\/19/ / 
and the Z)j's are given by, 

A = -9 [Xl,l(A Mi) + X3,l(A Mi) - 2X2,l(^,Mi) - 4(xe,5(^,w) + Xl6,5(^,W)) 
+8(x3,4(/3, Mi) + Xn,4(/3, Mi) + Xn,5(/3, Mi)) 

-16(X2,3(/3, Mi) + X8,3(/3, Mi) + X7,4(/3, //j)) 

+2 5 ( X5 , 3 (/3, Mi) + V2(X22jW, Pi) + V2X8,7(P, Mi))) 
-2 6 V2xi5, 7 (/3, Mi) 
-2 7 (xio,g(/3, Mi) + X28,g(/3, Mi)) 

+2 8 (xej(/3, Mi) + X5,7(/5, Mi) + Xi9,g(/3, Mi) + Xi9j(P, Pi) + X20j(P, Pi)) 
-2 8 (X8,nW, Pi) +Xio,n(/8, Mi) + X30,ii(A Mi) +X32,ii(/S,Mi)) 
+2 9 (x5,s(/3, Mi) + X2i,s(/3, Mi) + Xi9,n(0, Mi) + X2i,n(/3, Mi)) 
-2 9 (x3,5(/3, Mi) + Xi3,5(/3, Mi) + Xi2j(/3, Mi) + Xi3,7(£, Mi)) 
+2 10 ( X8i5 (/3, Mi) -Xi3, 8 (/3, Mi))] , (49) 
where [i = A,R). 

2.3.3 Calculation of £40 open diagram 

Performing the usual manipulations on the £40 term we obtain, 

Ao = \p 3 Jf[dR i f 2 (R l )v*(R 1 )h(R 2 )h(R 3 )h(R 4 )h(R 5 ) 
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Figure 7: Density dependence of the Em component. 



5(R 1 +R 3 + Ra)5(R 5 -R!- R 2 ) . 
The latter equation is simply, 

E40 = ~/J dR^iR^iR^vl^)} 2 , 
and after substitution of eq. ([1]) , 

vr 9 / V / c / 48 336 63 - 64^ 

40 ~ /3 6 \J3 \ 19^19 + 1697l3 + 294^7 + 
20^/6- 27 3 177 3 | 19jV3 \ \ 

150^5 22^n 256^ + 1024 + 576 ) + rVr aVa ) 

where, 

1 255 63 

p i = ^Xx,xiP,lH) ~ 25gX2,i(A/"i) + Y2gX3,l(AMi) 

+ ^X4,l(/ 3 > Mi) + |X3,2(/3,Mi) - ^X5,2(/3, /ij) 
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14 1 

8 14 

-^-/=X7,3(/3,Mi) + 2^3(/3,Mi) + g-^Xio,3(/3,MiJ 

1 

-4X11,3 (A Mi) - 8x7,6(/3, Mi) + 16xi3,e(/3, Mi) 
-8xi9,e(/3,Mi) 



5* 
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Pa(f™ 3 ) 



Figure 8: Density dependence of the £40 component. 



The density dependence of the E^q component is displayed in Fig. [51 



2.3.4 Calculation of E±c connected diagram 

In this subsection we estimate, 
lp 



1 p f 

Eac = ^ J n dr i v *( r ^)f 2 ( r ^) h ( r 23)h(r u )h(r 41 )h(r 2 4)h(ri 3 ) . 



(53) 
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Using the same technique as above we introduce 5 kernels, 

E 4C = -J I l[dR i v*(R 1 )f\R 1 )h(R 2 )h(R 3 )h(R 4 )h(R 5 )h(R 



6) 
1 

5(Re + R 3 + Ri)S(R 5 - R 2 - R 3 )S{R e - Ri - R 2 ) . (54) 

The identities 

(2vr) 3 5 {R 6 + R 3 + R 4 ) = J dqt exp(q 1 .( J R 6 + R 3 + R 4 )) , (55) 

(2tt) 3 J (R 5 -R 2 -R 3 ) = J dq 2 exp(q 2 .{R 5 -R 2 - R 3 )), (56) 

(2irf5 (R 6 -R 1 -R 2 ) = J dq 3 exp(q 3 .(Re - Ri - R 2 )) , (57) 

are introduced in Eq. (|54p . The latter equation is then expressed in terms of 
the Fourier transform of h, labelled h and that of v*f 2 , labelled G. More 
precisely, 

1 /• 3 ~ ~ ~ ~ ~ 

EiG = 2048tt 9 ^ J h (Ql ~ Q2)h(qi + q 3 )h{q 2 + q 3 )h(q 1 )h{q 2 )G(q 3 ) . 

i=i 

(58) 

For the sake of simplicity use is made of the scaling qj i— > 2y2/3qj in Eq. (|58p 
in order to eliminate the /3 dependence in h. Introducing, 

h*(q) = -2e- 2 * + ^~ q2 , (59) 

such that h(2\/2(3q) = (^/b) 3 h*(q), the equation ([55]) becomes, 

E 4 C = -3TT=3'° 3 / II d 9i ^(^l " Q2)h*(qi + q 3 )h*{q 2 + q 3 ) 

P b VK J i= i 

h*( qi )h*(q 2 )G(2V2Pq 3 ) . (60) 

The auxiliary integral 

r 2 

J (^3) = / II d( H h *(<H - Q2)h*( qi + q 3 )h*(q 2 + q 3 )h*( qi )h*(q 2 ) , (61) 
J i=l 

entering the definition of E 4c , has been performed by using the Cartesian 
coordinates. We obtain after calculations J(q 3 ) = vr 3 Q(q 3 ) with, 

9\/2 + 2V3 , „ 2 \ 8 ( 22 2 



0(9) = 72 exp(-2 q ) + —= exp ^ 
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32 2 \ / 1 1 \ / 3 2 
exp I — — q + = = exp 



19 v/19 V 19 J V56\/7 20 \/57 V 2 

2 / 10 2 \ 1 / 7 2 

exp ——q H t= exp —-q 



21 V21 \ 7 J 30 ^/30 \5 



I 2 \ / 4 2 \ 1 / 19 2 



1 " ^72tJ exp r 3 q ) + ^ exp r i6 q 

1 / 13 2 \ 9- 16 V3 2 



88 \/iT V U / 36864 



exp -— + ; exp(-g 2 ) . (62) 




Figure 9: Density dependence of the E^q component. 
Therefore, 



Eac 



4^0F 3 3 



P 



dqs Q(q3)G(2V2(3q 3 ) 



(63) 



is expressed in terms of the 3D Fourier transform (G) of v*(r)f 2 (r), which is 
given by, 



G(q) 



7T 



3/2 



V2 p 



exp 



r 



^3 eX P 



4 (5 2 



16 



I 



30 



20 



10 



-10 



-20 



-30 



- 

- 










E 4 
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Figure 10: £2, -E3, £4, E to t density dependence. 



+V R 



+ 



1 



4/4 

1 



Atf R (/3 2 + /4) 3/2 



(2 P+Vl )3/2 

2 



exp 



4(2 /3 2 + /4), 

-.2 \ 



exp - 



V A 



4(/3 2 + ^) 



exp 



4/4 



4/4 



(/3 2 + /4 ) 3/2 



exp 



+- 



1 



4(/3 2 + ^), 

/ „2 



(2/32 + ^)3/2 ex P^ 4 (2^ + Ai 2^ 
we are left with the final expression of E^c'- 



4V2 vr 9 / 2 3 
= — 

239063 

+ 



3142463 



200 y/2 657 y/2 
+ 



(3 \ 276480000 867 y/51 1225 



19 



+ 



243 



A+C r V r -CaVa 



■ (64) 



2074464 y/2 64 >/3 2450 >/35/ 
The latter equation is expressed in terms of the Cj's, i = A,R which are listed 
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Figure 11: Total energy for the case of the Ali-Bodmer potential compared to 
the Gogny Dl and DIN 

in the Appendix. 

The contribution to the energy supplied by each component of the cluster 
expansion is summarized in FigJIOIfor the AB interaction. Clearly the major 
contributions come from E2 and E3 whereas the four-body contributions are 
very weak. In contrast to two-body contributions the three-body and four- 
body contribution are repulsive and the equilibrium arises from a delicate 
balance between E2 and £3 components. 

We display in Fig. QT] a comparison of the EOS determined using the 
Ali-Bodmer potential with EOS obtained using Dl and DIN parametrizations 
of the Gogny potential. More details are given in Fig. [T21 While the equi- 
librium point predicted by the AB interaction is deep and lies at the normal 
nuclear matter density, Gogny parametrizations predict a shallower minimum 
at densities close to the Mott density. 

The four-body contribution to the total energy is small compared to other 
components. This is explained by an almost complete cancellation between 
E^r and E±d while the most connected diagrams E40 and E±c are intrinsically 
small due to increased number of bonds (see Fig. [T^|) . 
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Figure 12: E2,E^,E^,E dependence of the a-matter density p.(Ali-Bodmer 
and Gogny) 



3 Concluding remarks 

By assuming a simple but realistic functional form for the TBCFN's and em- 
ploying Gaussian- like potentials, we obtained analytical expressions for the 
first six terms (one corresponding to the two-body, one to the three-body 
and four for the four-body diagrams) occurring in the cluster expansion of an 
extended uniform system of structureless, indistinguishable a particles in a 
uniform background of neutralizing charge. The results are pointing to a sat- 
uration of the alpha matter at densities around the nuclear matter density (for 
Ali-Bodmer potential) and for the Dl interaction at almost one tenth of this 
density. We also demonstrated the strength of the folding method, previously 
applied in the calculation of the heavy-ion potential [12j to the calculation of 
cluster integrals of a many-boson system. 
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Figure 13: E^, E^, £40, E^c, E4 dependence of the a-matter density p.(Ali- 
Bodmer and Gogny) 



4 Appendix 

In this section we report the coefficients Cj entering the final expression of 
E±c- They are given below in terms of the coefficients, 



Xr 



(mf3 2 + ntff/ 2 



We have, 
Ci = 

+ 

+ 



129 1 1 
+ 7= + 



1 



2048 384 V3 24 

X3,2(/5,Mi) 



513 1 
4096 ~ 768 V3 ~ 12 ^6 

X3,l(£,W) 



(65) 



40 V5 7 Vn, 

+ (4^-7^71) X4 ' l(A ^ ) + (-5^ + iT7ii) X4 ' 3(/3 ' w) 

' V2 1 \ ,„ , ( 1 1 



+ 

+ 
+ 



J>/7 20 \/5, 
' \/2 1 



5^ 7 >/l4, 



X5,2(/3,^i) + 



40 \/5 7 v/14 



3 V6 



1 1 

+ 



5^ 14 ^ Xio,3(A^) + ^Xio, 7 (/3, W ) 

\/2 1 

+ Xi2,5(/3, I*) + 2 \/2 Xi3,n(/3, Mi) - ^~"7| Xl7 ' 5 ^' 
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- -^/q ^ 17 ' 7 ^' ~ ^ Xl9,16(/3, Mi) - g X22,13(/3, Mi) 

+ X24,7(/3, Mi) - 4 \/2 X24,ll(/ 3 ) Mi) + ^ X32,19(/3, Aij) 

+ 2 \/2 X35,ll(/S, m) + - X35,13(/3, Mi) + 2 X35,16(A Mi) 

- ^ X48,i3(/3, Mi) - v/2 X51,16(/5, Mi) 

- X5i,i9(/3, Mi) + ^= X70,ig(/3, //j). (66) 
where (i = A, R). 
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